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Social interactions vary in time and appear to be driven by intrinsic mechanisms, which in turn 
shape the emerging structure of the social network. Large-scale empirical observations of social 
interaction structure have become possible only recently, and modelling their dynamics is an actual 
challenge. Here we propose a temporal network model which builds on the framework of activity- 
driven time-varying networks with memory. The model also integrates key mechanisms that drive 
the formation of social ties - social reinforcement, focal closure and cyclic closure, which have 
been shown to give rise to community structure and the global connectedness of the network. We 
compare the proposed model with a real-world time-varying network of mobile phone communication 
and show that they share several characteristics from heterogeneous degrees and weights to rich 
community structure. Further, the strong and weak ties that emerge from the model follow similar 
weight-topology correlations as real-world social networks, including the role of weak ties. 


I. INTRODUCTION 

The emergent structure and dynamics of social in¬ 
teractions are consequences of individual and collective 
decision-making processes pQ. Understanding the driving 
forces behind social network formation has been a great 
challenge that has attracted lots of attention during the 
last decades. Even though conventional experimental ap¬ 
proaches (based on e.g. surveys) have revealed fundamen¬ 
tal rules of social behaviour [2], it has not been possible 
scale such approaches up to sizes required for observing 
features of large-scale social structures, let alone their de¬ 
tailed dynamics. However, technological advances have 
helped to overcome these limitations through detailed 
digital recording of social interactions ehh- Following 
these advancements, a lot of effort has been put on map¬ 
ping social structures from digital traces, and on estab¬ 
lishing methodology for the analysis of social networks 
connecting millions of individuals. These efforts have led 
to discover the small-world architecture [5] , the heteroge¬ 
neous network structures induced by preferential attach¬ 
ment mechanisms mm, the different roles of strong and 
weak ties 0! and their relationship to triangle forma¬ 
tion m and the emergence of communities mm- 

The earliest data-driven studies of social networks were 
mainly focused on static network structure. However, it 
is evident that this approach both unnecessarily discards 
available information - electronic interaction traces typ¬ 
ically come with time stamps - as well as hinders de¬ 
tailed understanding of the mechanisms shaping network 
structure. In the recent years, steps have been taken to 
come around these limitations, with increased focus on 
both network dynamics (changes in link structure, such 
as in )13| ) and temporal network structure [14] (time- 
stamped interaction events between individuals giving 
rise to e.g. temporal motifs BSD- In particular, stud- 
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ies of temporal networks of human interaction have re¬ 
vealed tie dynamics with memory and reinforcement pro¬ 
cesses mm, bursty dynamical patterns PM9] , tem¬ 
poral homophily in communication motifs m as well 
as different strategies of communication and tie mainte¬ 
nance pT] , 

Modelling temporal networks and their co-evolution 
with processes unfolding on the networks is still one of the 
key challenges in the field. Empirical studies employing 
null models with shuffled interaction sequences |??rf2T| 
have pointed out that bursty interactions together with 
features such as weight-topology correlations strongly in¬ 
fluence the speed of any type of information diffusing on 
the temporal network El HIM!. However, a limiting 
factor of this approach is that it only allows to remove 
existing temporal inhomogeneities and correlations from 
empirical sequences in order to understand their effects. 
Because of this, the null model approach does not pro¬ 
vide ways to continuously adjust the level of inhomo¬ 
geneities or temporal correlations. This limitation can 
be overcome by resorting to artificial models of tempo¬ 
ral networks |28] [25] that allow to fine-tune the level of 
different features one by one. One promising direction is 
the activity-driven framework [29] that assumes heteroge¬ 
neous activity of agents. It provides a general foundation, 
which can be extended by incorporating additional tem¬ 
poral and structural correlations. Studies of the activity- 
driven time-varying network model have demonstrated 
that memory effects can explain the emergence of strong 
and weak ties in social networks m, while endo- and 
exogeneous tie formation processes may control the evo¬ 
lution of business alliance networks [30]. The same model 
has been used to study how time-varying interactions in¬ 
fluence the dynamics of random walkers ED EH , rumour 
spreading m, and epidemic processes pi3]. 

In this paper, our aim is to consider further social 
mechanisms that would make models of time-varying so¬ 
cial networks more realistic. We introduce a temporal 
network model with adjustable community structure and 
emergent weight-topological correlations via the exten- 
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sion of the activity-driven time-varying network model 
with three additional mechanisms. These are i) rein¬ 
forcement processes to model memory-driven interaction 
dynamics of individuals m, ii) focal and cyclic closure 
motivated by the work of Kumpula et. al. [T2] [31, to 
capture patterns responsible for the emerging commu¬ 
nity structure, and iii) a node removal process. Using 
this temporal network model we demonstrate the effect 
of the scalable community structure and social reinforce¬ 
ment on information spreading, which co-evolves with 
the time-varying interactions. After the introduction of 
the model, we discuss its temporal behaviour, and match 
the emerging network properties with empirical observa¬ 
tions. We discuss the present higher order correlations 
and their effects on spreading processes, and finally we 
conclude our work and discuss possible future directions. 


II. MODEL 

Our model definition builds on the activity-driven 
time-varying network (ADN) model introduced by Perra 
et al. 29]. The original model takes N nodes that are 
each assigned a priori with an activity probability per 
unit time ai = ijXi. Here Xi is the activity poten¬ 
tial of node i, drawn from a suitable distribution F(x{) 
{xi £ [e, 1 ] bounded by a minimum activity e), and rj is 
a time rescaling factor determining the average number 
of active nodes per unit time to be rj(x)N. At the be¬ 
ginning of each time step (time t), the model network 
Gt = ( V t ,E t ) is initially disconnected ( E t = 0) and it 
evolves via the following steps: (i) Through random se¬ 
quential update each node i becomes active with prob¬ 
ability aiAt and connects to to other randomly selected 
nodes in G t . (ii) At the end of each step (time t + At), 
all created links are removed and the process moves on 
to the next iteration. Note that during a single step, a 
node can participate in temporal interactions by creat¬ 
ing a link (if it is active), or by receiving link(s) from 
other active nodes. The activity potentials of humans 
are commonly observed to be heterogeneous [55U3T] . We 
approximate this with a power-law activity potential dis¬ 
tribution F{xi) ~ x with exponent 7 = 2.8, based on 
empirical observations and earlier modelling work m 
Furthermore, since we aim to model one-to-one interac¬ 
tions, typical for mobile telephone communication, we set 
to = 1. In addition, without loss of generality, we fix the 
parameters 77 = 1, e = 10 -3 and At = 1 [171 . 

A. Social mechanisms 

The first mechanism to be integrated with the activity- 
driven model is memory in terms of frequently repeated 
interactions between a node and its peers who have been 
contacted before. This has been considered earlier by 
Karsai et al. El with reinforcement processes, where a 
node remembers its connected neighbours, and depend¬ 


ing on its degree n at f, either creates a new link with 
probability p(n) = c/{n + c), or interacts with a neigh¬ 
bour with probability 1 — p(n) = n/(n + c) to reinforce 
an existing tie. The introduction of memory defines a 
non-Markovian dynamics where past interactions influ¬ 
ence the actual decision of a node. The larger a node’s 
egocentric network size n, the smaller the chance for it to 
create a new link in the actual iteration. This captures 
the realistic assumption that interactions in a social net¬ 
work are not random but favoured with actual friends 
whose number is strongly limited by cognitive capaci¬ 
ties [3H132] ■ This decision-making mechanism naturally 
leads to the emergence of weak and strong ties, weight 
heterogeneities, and a more realistic degree distribution 
in the integrated network structure El- Note that in 
this formulation c scales the strength of memory, which 
may capture the attitude being a social keeper or social 
explorer [211 . Here, for simplicity, we set c = 1 for each 
node and leave the exploration of the effect of varying 
memory strengths for future studies gD]. 

The second mechanism we consider is that of tie cre¬ 
ation with closure processes. Here, cyclic closure, the 
tendency to form network triangles, shapes the social 
structure at mesoscopic scale, and is responsible for the 
emergence of communities m- The mechanism of fo¬ 
cal closure, on the other hand, is independent of network 
structure and represents the formation of ties between in¬ 
dividuals with shared attributes or interests. It is driven 
by the propensity to seek cognitive balance between con¬ 
nected egos m Eg as suggested by earlier theories in so¬ 
ciology mm- An applicative definition of cyclic and 
focal closure in general is given by Kumpula et al. mm , 
who model cyclic closure as biased local search. To the 
contrary, focal closure is modelled as an unbiased global 
random search. 

Finally, the third ingredient of the model is the process 
of node removal. It allows the model network to reach 
an equilibrium state where its overall structural charac¬ 
teristics become invariant of time. 


B. Model definition 

Here we introduce a time-varying network model, 
which integrates memory and reinforcement processes 
with closure mechanisms in the activity-driven model def¬ 
inition. Our ultimate goal is to provide a temporal net¬ 
work model with adjustable community structure and 
weight-topology correlations to understand their role in 
shaping the emergent network structure and co-evolving 
dynamical processes. 

Our model takes N initially disconnected nodes and 
updates them in random order at each iteration step. 
In one iteration step, a node may become active with 
probability ai = pxi or be deleted with probability pd- 
Note that to keep the network size constant, for each 
deleted node, we add a new node to the network in the 
next iteration step. While a node i is alive, it remembers 
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its already connected neighbours j £ Vf and the weight 
v/-,: of interactions with each of them. 

L J 

If the node i becomes active, then, depending on its 
actual degree i%i, it can attempt to form a new link with 
probability p(rii) = 1/ (n, + 1), or otherwise interact via 
an existing link with probability (1 — p(ni)). In the lat¬ 
ter case it randomly selects one of its neighbours j with 
probability pfj = w\j/J^kevj w lk weighted by the num¬ 
ber of their interactions. The two nodes then interact 
and increase their link weight w\j by 6 (reinforcement 
process). 

On the other hand, if the node decides to form a new 
link, it may follow different scenarios. In all cases, the 
new tie will initially have unit weight w* = 1. If the 
degree of the focal node is 0, it randomly picks another 
node from the entire network j (focal closure ) and forms 
a tie. Otherwise, it attempts to create a new link with 
the triadic closure mechanism. First, it chooses one of 
its neighbours j randomly with a weighted probability 
pfy If j has no other neighbours than i, node i looks for 
another random node to interact with ( focal closure) and 
forms a link. Otherwise, it looks for a random neighbour 
k of j (i ^ k) with a weighted probability pj k . If k is 
not an already existing neighbour of i (k ^ Vf) , the two 
nodes interact with probability pa, and close the triad 
by forming a link. Otherwise, with probability 1 — pa, 
node i follows the focal closure mechanism and instead 
forms a link with a randomly selected node (other than 
j and k). Finally, if k is already a neighbour of i, that is, 
k £ Vf , the two nodes interact and increase the weight 
of their existing link by S ( reinforcement process). At 
the end of each iteration step, all nodes finish their ac¬ 
tive interactions but remember their weighted egocentric 
network. For a pseudocode version of the algorithm, see 
Appendix A. 

III. MODEL NETWORK ANALYSIS 

In addition to the activity-driven model parameters 
whose values are fixed, our model has three intensive pa¬ 
rameters, pa, Pd, and S. By varying these parameters, 
one can simulate a rich variety of time-varying networks 
with several emergent structural properties and correla¬ 
tions. In the following, we explore how the properties of 
the emerging network structure depend on time and on 
the intensive parameters, and whether those properties 
match with those of an empirical temporal network of 
mobile phone communications. The dataset we use here 
contains 633,986,311 time-stamped mobile call commu¬ 
nication records of 6, 243, 322 customers of a mobile oper¬ 
ator (market share « 20%) in a European country. Cus¬ 
tomers are represented as network nodes, connected via 
16,783,865 weighted mutual links with weights defined 
as the number of calls between customers who mutually 
called each other at least once during the observation 
period. The presented network properties of the mobile 
phone call (MPC) network were calculated from a static, 


aggregated representation of the social network structure 
obtained by integrating the temporal interactions over 6 
months. 

In the following, model networks were generated via 
large-scale numerical simulations with N = 10,000 
nodes, and results were averaged over 100 independent 
realisations (if not noted otherwise). For each realisa¬ 
tion, we measured the network parameters by consider¬ 
ing links that are actually present in the network, i. e. we 
disregarded links of removed nodes. 

A. Temporal features 

The introduced network model is inherently tempo¬ 
ral and simulates time-varying interactions between in¬ 
dividuals. To explore its overall temporal behaviour, we 
measured two general network properties as a function of 
time. The first property is the average degree defined as 
(k)(t) = N~ x ki(t), where the sum runs over all nodes 
i, and ki denotes the number of established connections 
of node i at time t (i.e. ki(t) = \Vf\). Second, we mea¬ 
sure the average local clustering coefficient C[t), defined 
as the fraction of the real and possible numbers of trian¬ 
gles around a node given its number of links, averaged 
over the whole network [44 j . C quantifies the density of 
triangles in a network; in social networks, the existence 
of communities typically gives rise to high triangle den¬ 
sities. 





FIG. 1: Average network properties as the function of time. 
We depict the evolution of a) the average degree ( k)(t ), b) 
clustering coefficient C, and c,d) the fraction of events 
creating ( n cr (t )) or reinforcing (n r /(t)) a link with various 
values of the deletion probability pd (for exact values, see 
legend). For each calculation, we fixed pa = 0.5 and <5 = 1. 

As follows from the model definition, the process starts 
from a set of disconnected agents, thus the measured 
properties are trivially zero at time t = 0. However, 
as time goes by and ties are formed via temporal interac¬ 
tions, ( k)(t) starts increasing until the network reaches a 
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FIG. 2: Demonstration of the emerging structure in the time-varying network model. Panels (a-c) depict simulated networks 
with fixed 5 = 1 and varying pa = 0.5, 0.9, and 0.995 (pd = 4e — 5, 2e — 5, and 1.04e — 5) respectively. Panels (d-f) shows 
networks with fixed pa = 0.995 with varying 5 = 0, 0.5, and 1.5 (pd = 3.5e — 5, 1.7e — 5, and 8e — 6) respectively. Each panel 
depicts the actual structure of a network with N = 500, in its stationary state. Links are coloured according to their weight 

(darker link colour = stronger link weight). 


stationary state with constant average degree. The time 
it takes to reach this equilibrium state strongly depends 
on the choice of the node deletion probability pd as shown 
in Figjlja. If pd is too small, nodes remain in the sys¬ 
tem for a long time and even nodes with small activity 
levels have time to evolve their egocentric network re¬ 
sulting in a slow relaxation to the stationary state. On 
the other hand, as pd is increased due to the finite life¬ 
time of nodes, less active nodes are removed before their 
egocentric structures are fully evolved. Because of this, 
the network reaches equilibrium faster. This explains the 
decrease of stationary average degree with increasing pd- 
The average degree can be tuned approximate the em¬ 
pirical average degree of the MPC network (dashed line 
in Fig{l]a) 

Measurement of the time dependence of the the av¬ 
erage clustering coefficient (see Figjljb) yields a fairly 
similar scenario. However, here, after the initial increase 
of C{t) it reaches a maximum, followed by a decrease 
and relaxation to a stationary value. This is because 
triangles tend to emerge right in the beginning of the 
process, followed by the evolution of strong ties. After 
nodes have created their first links and closed triangles 


between them, interactions begin to frequently take place 
on existing links that are reinforced in the process. Thus 
the local search of active nodes is biased towards strong 
ties that are part of already formed triangles. In late 
times nodes attempt less frequently to create new trian¬ 
gles through the weak links that emerge throughout the 
dynamics by focal closure. Again, there is saturation due 
to node removal influenced by the rarely active but sur¬ 
viving nodes, who keep introducing random links in the 
network. By choosing an appropriate parameter values, 
the clustering coefficient can again be tuned to compa¬ 
rable to its corresponding empirical value (dashed line in 
Fig{T}b). 

The above picture is fully supported by measurements 
capturing the fraction of events which create new links 
n cr (t ) = #e cr ./(#e cr .(t) + #e r /(t)), or reinforce already 
existing ones n r f(t) = #e r //(#e cr (t) + #e r /(t)). Here 
#e cr (#e/ r ) denotes the number of such events at time t. 
As seen in FigjTJc these measures take values of 1 and 0 
(respectively) at time t = 0, as all events initially create 
new links. However they rapidly approach a constant 
value as it is visible in FigjTJd. On the other hand, as 
Pd is increased, the number of events reinforcing already 
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existing links is increased, while events responsible for 
the creation of new links become less frequent. This is in 
accordance with our interpretation above. The measures 
can be tuned to reflect empirical network values here as 
well (black lines). 

In the following, if not noted otherwise, we set pd = 
0.5 x 10 -5 and run the simulations for t = 150, 000 time 
steps to ensure the model networks to always reach their 
stationary state with average properties invariant in time. 


B. Static features 

Next we concentrate on the effect of different mech¬ 
anisms on the emerging network structure in the sta¬ 
tionary state. In Fig(2j we visually demonstrate the role 
of the cyclic closure (a-c) and link weight reinforcement 
(d-f) mechanisms in the emerging model networks. In 
panels Figj2](a-c), we have kept the reinforcement incre¬ 
ment constant, <5=1, while varying the cyclic closure 
probability pa between 0.5 and 0.995. The node dele¬ 
tion probability pd has been chosen to yield networks 
with suitable link density for visualisation (for exact pa¬ 
rameters see the figure caption). When pa is small, the 
emerging network structure is densely connected and ap¬ 
pears more like a random structure since link creation 
is driven by focal closure. Nevertheless, because of the 
reinforcement process, weight heterogeneities emerge al¬ 
ready in this case. More interestingly, by increasing pa, 
communities are seen to emerge, with weight-topology 
correlations that are in line with the Granovetterian pic¬ 
ture Bang, where strong ties connect nodes inside com¬ 
munities, while weak ties emerge between them (darker 
link colour = stronger link weight). However, focal clo¬ 
sure alone is not sufficient for the emergence of commu¬ 
nity structure as shown in panels Fig{2](d-f). Here, the 
triangle formation probability has been kept constant, 
Pa = 0.995, while the reinforcement parameter 5 has 
been varied between 0 and 1.5. Even though a triangle 
would almost always be closed by the local search if a 
suitable node were found, without link weight reinforce¬ 
ment (d), the local search will be disperse and often re¬ 
sulting in focal closure because hitting a neighbour with 
no other neighbours. However, by increasing S the local 
search will get more biased towards emerging strong ties 
and lead to the evolution of tight community structure. 
Note that since S controls the strength of the local search 
mechanism, it also scales the size of the communities. 

To quantitatively explore the emerging network struc¬ 
ture as the function of the intensive parameters, we have 
measured the behaviour of general network properties. 
Fig|3] shows how the average degree of the model net¬ 
works depends on pa and <5 (panels a and b, respectively) 
for different deletion probabilities pd■ The average de¬ 
gree of the empirical MPC network is displayed for com¬ 
parison (dashed horizontal line). In each case, the av¬ 
erage degree decreases when the parameters controlling 
the probability of cyclic closure pa and the amount of 


link reinforcement 5 are increased, indicating that these 
parameters together with pd strongly control (k). How¬ 
ever, while the dependence of (k) onpA is concave, ( k)(5 ) 
is a convex function for any pd . Moving beyond aver¬ 
ages, panels c) and d) display the behaviour of the degree 
distributions P(k) for various pa and <5, indicating that 
there are strong degree heterogeneities. While the av¬ 
erage degrees strongly depend on these parameters, the 
overall distribution shapes do not change much. The dis¬ 
tribution of the empirical MPC is again displayed for ref¬ 
erence, showing a very similar shape than those produced 
by the model. 




FIG. 3: General network measures as a function of the 
intensive parameters, (a) and (b) show the dependence of 
the average degree (k) on pa and 5, respectively (for the 
former, 5=1 and the latter, pa = 1). Various values of pd 
have been used (see legend in (b)). The dashed black line 
depicts the average degree (k) of the empirical MPC 
network, (c) and (d) depicts the degree distributions P(k) of 
model networks with varying pa and 5, while (e) and (f) 
show the corresponding weight distributions P(w) . In 
panels (c) and (e) ((d) and (f)), we kept 5 = 1 (pa = 1) and 
Pd = 5 x 10 -5 . Black circles denote the corresponding MPC 
distributions. 

The strength of social ties, measured here as the num¬ 
ber of dyadic interactions, is commonly observed to be 
distributed heterogeneously in social networks. This is 
also the case with our empirical MPC network (black 
circles in Fig{3]e and f). Tie strength heterogeneity is an 
emergent property of our model networks where it comes 
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from the preference of nodes to reinforce existing links. 
This mechanism is independent from the search strate¬ 
gies, as evident in Fig(3je, where the model weight distri¬ 
butions are invariant to the selection of the cyclic closure 
probability pa (S and pd are kept constant here). The 
emerging weight heterogeneities are the consequence of 
the intrinsic memory process, which inclines nodes to re¬ 
inforce existing links rather than create new ones. There¬ 
fore, the weight heterogeneity naturally depends on the 
reinforcement parameter <5 (see Fig[3]f), which in turn 
mainly scales the tail of P(w) but does not affect not its 
functional shape, which appears to be a power-law with 
a constant exponent. 


C. Higher-order correlations 

Other than the observed degree and weight hetero¬ 
geneities in the model networks, certain parameter ranges 
result in the emergence of rich community structure. In 
real social networks, communities can be characterised by 
higher-order correlations as they are built of closed trian¬ 
gles with unevenly distributed links of various strength in 
the structure. Stronger ties that are maintained by fre¬ 
quent interactions tend to connect nodes inside communi¬ 
ties and shape the structure locally, while weak links with 
infrequent interactions are situated between communi¬ 
ties and keep the social structure globally connected. To 
show whether the modelled temporal networks can re¬ 
cover these characteristics, we have performed three sets 
of measures and compared our findings to the empirical 
results. 

As we have seen earlier, both the values of pa and <5 
influence the community structure. However, we expect 
that pa plays a dominant role here as it controls the tri¬ 
angular closing mechanism in the network. If pa = 0, 
links are created randomly and the clustering coefficient 
is very small as seen in Fig|4ja, where a constant (5 = 1 
were set for each measurement. By increasing pa , cyclic 
closure becomes more dominant, reflected in an increas¬ 
ing C with very high values as pa goes to its maximum. 
The clustering coefficient of the MPC network is depicted 
as the horizontal dashed line. There are several sets of pa¬ 
rameter values for which the model yields networks with 
similar clustering coefficient values, e.g. pa — 0.5, S ~ 1, 
and pd — 5 x 10~ 5 (see the crossing points of the black 
dashed and green solid lines in Figj4ja). Not surprisingly, 
5 plays a weaker role in the emergence of triangles. The 
reinforcement mechanism introduces a bias in the local 
search for new ties, which leads to tighter communities 
and more closed triangles as we increase 6, reflected by 
the slightly increasing C in Figj4]b. 

Second, we have measured weight-topology correla¬ 
tions to check if the model networks are constructed ac¬ 
cording to the Granovetterian weak-tie structure. In his 
seminal paper [8j Granovetter suggests that the fraction 
of common friends of connected individuals is positively 
correlated with their tie strength (i.e. their link weight 
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FIG. 4: Measures of higher-order correlations in the 
empirical and model networks. Panels (a) and (b) show the 
dependence of the average local clustering coefficient C on 
the parameters pa and 5, respectively, for the empirical and 
model networks. Weight-topology correlations are shown as 
the function of the same parameters ((c) and (d), 
respectively) by measuring the average overlap (O) as the 
function of cumulative tie strength P S um(w). Panels (e) and 
(f) depict the average local clustering coefficient of the 
residual networks after removing a fraction fo of nodes 
ranked by their overlap. Empirical results are shown in each 
panel with black dashed lines or circles, while measures on 
model networks are depicted by coloured solid lines with the 
corresponding parameters in the legend. In each 
measurements in panels (a,c,e) we kept 5 = 1.0, in panels 
(b,d,f) pa = 0.5, while in panels (c-f) pd = 5 x 10 -5 . 


here). The fraction of common friends can be measured 
by the link overlap [5] defined as Otj = riij/((ki — 1) + 
(kj — 1 ) — riij) where is the number of common neigh¬ 
bours of nodes i and j, and ki and k 3 are their degrees. 
To quantify weight-topology correlations, we measure 
the average link overlap (O) as a function of cumula¬ 
tive tie strength P CU m{w) that measures the fraction of 
links with tie strength smaller than w jjjj. In the MPC 
network, this function reflects positive correlations be¬ 
tween overlap and tie strength (black circles in Figj4jc 
and d) in accordance with earlier observations [9]. If the 
model networks have their triangle-formation mechanism 
or weight reinforcement turned off (pa = 0 or S = 0), 
no such correlation is found. However, for any positive 
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values of pa and S, a positive weight-topology correla¬ 
tion emerges in accordance with the empirical observa¬ 
tion and the hypothesis of Granovetter. For larger pa 
(with constant <5 = 1 ), the function is shifted as more tri¬ 
angles evolve, which is increasing the average overlap (see 
Figj4]c), while by increasing 6 (and constant pa = 0.5) 
(see Figj4jd) the correlation becomes stronger as ties with 
larger strength appear in the networks. 

At last, we have checked how the clustering coefficient 
C changes by removing a fraction of links fo in increas¬ 
ing order of overlap. The resulting functionf C(fo) can 
be divided into two regimes. Trivially the removal of 
links with 0 = 0 , which connect communities, does not 
decrease the number of triangles in the network but only 
decreases the degrees of nodes by removing links not par¬ 
ticipating in triangles. Consequently, by the removal of 
these links the clustering coefficient must increase. On 
the other hand, the removal of links with O > 0 that 
connect nodes inside communities reduces the number of 
triangles, resulting a decrease in C. This behaviour is 
observed in the empirical system (black circles in Figj4je 
and f) and also for the model networks. In Fig@ e as 
Pa 0 no triangles evolve in the network, and thus 
most of the links have 0 = 0 , and the first regime is 
extended. In the other extreme, as pa —> 1, most of the 
links have non-zero overlap and the second regime domi¬ 
nates. The best match with the empirical curve appears 
when pa — 0.5 (red curve in Figj4je) and ~ 40% of links 
evolve with 0 = 0. Note that the second regime of this 
function is very sensitive to the fine-grained structure of 
the actual communities, which causes the difference be¬ 
tween the model and empirical curves. In addition, by 
keeping pa = 0.5, we have checked the 8 dependence of 
this function (see Figj4jf). As the tie reinforcement 5 in¬ 
creases, links form tighter communities and the clustering 
coefficient increases in this case. 


IV. EFFECTS ON SPREADING PROCESSES 

We conclude our modelling study by demonstrating its 
capacity in simulating dynamical processes co-evolving 
with the temporal network structure with scalable in¬ 
tegrated social mechanisms. For this reason, we have 
selected the simplest possible spreading process, the 
susceptible-infected (SI) model [35], where each node can 
be in one of two mutually exclusive states; susceptible (S) 
or infected (/). Initially, every node is in the state S, ex¬ 
cept for a random seed with state /. The infection passes 
from any node in the / state to any node in the S state 
via a temporal interaction, but independent of the direc¬ 
tion of the actual temporal link [22] . In order to assure 
that the network is in the stationary state we initiated 
the spreading process after the temporal network of size 
N = 10,000 has evolved for t = 50,000 iterations. Af¬ 
ter this time we measured the i(t) = I(t)/N fraction of 
infected nodes until t = 150,000, where the process has 
reached maximal penetration (see Fig[5|a and b). Note 


that maximal penetration is not converging to 1 as there 
are always new susceptible nodes introduced in the net¬ 
work, moreover the temporal structure is not necessarily 
connected. For simplicity, we set the probability of trans¬ 
mission per interaction event to A = 1 (note that then the 
speed of transmission between two nodes is only limited 
by the frequency of their interactions). Further, we con¬ 
sider interactions as undirected and therefore they may 
transmit the interaction both ways. 



FIG. 5: Susceptible-Infected process co-evolving with the 
time-varying network model (a, b) Fraction of infected nodes 
as a function of time for different values of pa and 8 (for 
exact parameters see legend), (c, d) Time t$ o% of 50% 
penetration for different values of pa and 8. While in case of 
panels (a, c) we have kept <5 = 1, for panels (b, d) we have 
used pa = 0.99. Each simulation was run for t = 150, 000 
iteration steps with p d = t x 10~ 5 on networks of size 
N = 10,000. Results were averaged over 100 realisations, 
where we initiated the process from a single random seed at 
t = 50, 000. 


The simulation results depicted in Fig[5]a and b show 
that both the spreading process becomes slows when ei¬ 
ther the cyclic closure probability pa or the weight rein¬ 
forcement parameter 8 are increased. While in the first 
case this is mainly due to the emerging tight commu¬ 
nities, in the second case weight-topology correlations 
also play a role. The same conclusions can be drawn 
by measuring the speed of spreading, quantified as the 
time £ 50 % when the process has reached 50% penetra¬ 
tion. Tight communities evolve as pa —> 1, which con¬ 
strains the spreading process and explains the convex 
shape of the curve in Fig[5]c. On the other hand, weight- 
topology correlations appear even for small values of <5 
(as we have shown in Fig@ c). These immediately slow 
down the speed of spreading and cause the concave shape 
of the function in Fig[5]d. Thus, one can conclude that 
even though more pronounced community structure has 
considerable effects on the spreading process, already 
weak weight-topology correlations can strongly affects 









the speed of spreading. 


V. CONCLUSIONS 

The temporal network representation takes into ac¬ 
count the time-varying nature of interactions between 
entities instead of considering only the static network 
structure of their connections M- So far, this approach 
has mainly focussed on analyzing empirical data, while 
there is still a lack of appropriate network models. In this 
paper we introduce a model of temporal social networks 
with emerging heterogenous structure, communities, and 
higher-order correlations. As a starting point, we take 
the recently introduced modelling framework of activity- 
driven networks with heterogeneously distributed activ¬ 
ities of individuals [29] and memory processes |l7j . We 
extend this model with three mechanisms that have been 
earlier used for generating static weighted networks m, 
such as social reinforcement , cyclic and focal closure. 
These mechanisms arguably drive the tie formation of 
individuals and are responsible for the emergence of com¬ 
munities and the global connectedness of the network. In 
addition, we introduce a node removal process allowing 
the model network to reach an stationary state, where 
its overall structural characteristics become invariant of 
time. By adjusting the model parameters, it is possible 
to recover several characteristics of a real-world tempo¬ 
ral network of mobile phone calls. These include hetero¬ 
geneous degrees and weights, rich community structure, 
and weight-topology correlations and weak tie properties 
as suggested by Granovetter [8]. 

The main advantages of this model are that (a) it is 
able to mimic the time-varying nature of interactions 
characterising several real systems; (b) it helps to under¬ 
stand the importance of different mechanisms in shap¬ 
ing the emerging network structure; and (c) it allows us 
to test their effect in a scalable fashion on co-evolving 
dynamical processes. This last point is of special im¬ 
portance, as e.g. spreading processes have been shown 
to have significantly different critical behaviour when co¬ 
evolving with time varying interactions. In addition, by 
varying the model parameters one can control the average 
degree, strength of ties, size of communities, clustering, 
and interconnectedness of the network. These features 
can control e.g. the speed of spreading as we demon¬ 
strated via a simple study on SI processes, but play cru¬ 
cial effects in case of other processes like random walks, 
the diffusion of information, epidemics, or social conta¬ 
gion processes, making the model an ideal testbed for 
simulating such processes. 
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Appendix A: Algorithmic model definition 

Here we describe the algorithmic implementation of 
the temporal network model introduced in the main text. 
The main function Temporal network model(Gt=o, 
T, pA; Pd, 5) takes as input a weighted network with 
empty link set G t= o = (V, 0, {at}, {w^}) with size N, 
pre-assigned node activity probabilities per unit time 
{ai}, and initiallinkweights{wlj = 0} . Further pa¬ 
rameters are T: the number of iterations; pa- proba¬ 
bility of triadic closure; Pd : probability of node deletion; 
and S: link reinforcement increment. In the pseudo code, 
Ui = Vf denotes the current degree of node i; randif) is 
a pseudo-random number generator, which returns a ra¬ 
tional number between [0,1]; and = w\^/ w ij 
is the probability to select a random neighbour j from 
the current neighbour set Vf of node i weighted by the 
number of interactions between them performed since the 
link was created up to time t (ffZjev* Pfj = !)■ Temporal 
events between nodes i and j performed at the current 
iteration step t are denoted by [t,i,j]. 


ALGORITHM 1: Temporal network model 


Input : G t = o , T, pa, p d , S 

for t = 1 to T do 

for N randomly selected node i £ V do 
if rand () < p d then 

remove all links of node i o delete and re-insert i 

end if 

if rand( ) < a; then > node i is active 

if rand() < pfnf) then o link creation 

if m = 0 then > node i has no neighbour 

Focal closure(i, {i} ■ Gt) 
else 

Cyclic closure(i, Gt) 

end if 

else > link reinforcement 

Select a j neighbour of i with probability pfj 
Event [t,i,j] 
w\j T — d 

end if 
end if 
end for 
end for 


The main function calls two subroutines called Cyclic 
closurefi, Gt) and Focal closurefi, Gt, X) with input 
parameters i denoting the currently active node; Gt the 
current weighted network structure; and X exception set 
of nodes to connect in the current call. 
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ALGORITHM 3: Focal closure : as input i is the current 
active node; X is the exception set of nodes to connect; and 
Gt is the current weighted network 


Input : i, X, G t 

Select a random node j € V \ X 
Event [t,i,j] 
w\j = 1 


if rand() < pa then 
Event [t,i,k\ 

Wik 1 

else 

Focal closurefj, {i,j, k} LiVt, Gt.) 

end if 
else 

Event [t,i,k\ 

W ik~^~ = $ 

end if 
end if 


ALGORITHM 2: Cyclic closure: as input i is the current 
active node; and Gt is the current weighted network 


Input : i, Gt 

Select a j neighbour of i with probability pfj 
if rtj = 1 then > node j has one neighbour 

Focal closurefj, {i,j} U V) , Gt) 
else 

Select a k(^ i) neighbour of j with probability pj k 
if (i,k) £ Et then > triangle ( i,j,k ) is not closed 
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